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ABSTRACT 

In the spirit of classic works of Wilson on the renormalization group and operator product expan- 
sion, a new framework for the study of the theory space of euclidean quantum field theories has been 
introduced. This formalism is particularly useful for elucidating the structure of the short-distance 
expansions of the n-point functions of a renormalizable quantum field theory near a non-trivial fixed 
point. We review and apply this formalism in the study of the scaling limit of the two dimensional 
massive Ising model. Renormalization group analysis and operator product expansions determine all 
the non-analytic mass dependence of the short-distance expansion of the correlation functions. An 
extension of the first order variational formula to higher orders provides a manifestly finite scheme for 
the perturbative calculation of the operator product coefficients to any order in parameters. A pertur- 
bative expansion of the correlation functions follows. We implement this scheme for a systematic study 
of correlation functions involving two spin operators. We show how the necessary non-trivial integrals 
can be calculated. As two concrete examples we explicitly calculate the short-distance expansion of 
the spin-spin correlation function to third order and the spin-spin-energy density correlation function 
to first order in the mass. We also discuss the applicability of our results to perturbations near other 
non-trivial fixed points corresponding to other unitary minimal models. 
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1 Introduction 



The equivalence of statistical mechanics models and quantum euclidean field theories with an ul- 
traviolet cutoff has provided us with the exchange and deeper understanding of very fruitful notions, 
such as the renormalization group (RG) and spontaneous symmetry breaking, between the two disci- 
plines. The modern renormalization group language offers a beautiful qualitative description of how 
the existence of second order (continuous) phase transitions in a statistical mechanical system is rele- 
vant in defining a local quantum field theory In particular the existence of a critical point and the 
number of parameters one needs to fine-tune in order to reach the critical point of a statistical model 
correspond to the existence of an ultraviolet fixed point and the number of relevant operators of the 
continuum field theory, respectively. In fact the exact solvability of statistical lattice models gives us 
detailed quantitative information on the existence and proper definition of the continuum limit of a 
quantum field theory and, in principle, allows for a non-perturbative evaluation of all the correlation 
functions in the theoryQ 

In many situations where the lattice model is not exactly solvable, the exact solvability of the 
fixed point theory is expected to allow for perturbative calculation of the correlation functions of 
the model near criticality. This is the situation ordinarily encountered in field theory calculations 
where the fixed point theory is trivial, i.e., a free field theory. It is reasonable to ask whether the 
same is true of non-trivial exactly solvable fixed point theories. Fixed point theories are the fixed 
points of the renormalization group action on the theory space of all quantum field theories. Due to 
the complexity of theory space with respect to the RG flows, information about theories not in the 
immediate neighborhood of a given ultraviolet fixed point is accessible only non-perturbatively. It is 
believed that, in many physically interesting cases, under the action of the renormalization group the 
theory flows from an ultraviolet (UV) fixed point to an infrared (IR) fixed point in the theory space. 

It is very difficult to give a precise definition of the space of all quantum field theories. In fact this 
definition is dependent on the field-theoretic issues under consideration. For a better understanding 
of theory space it is clearly important to have a perturbative scheme for studying off-critical theories, 
in the vicinity of and away from their non-trivial fixed points, given the information encoded in the 
fixed point theories only. When off-critical theories are defined as perturbations of critical theories 
we should in general consider perturbations by relevant, marginal, and also irrelevant operators of the 
critical theory. Among these, the perturbations which are integrable are of special interest. Comparison 
of perturbative calculations in such models with exact results may shed light on how non-perturbative 
information can be extracted from such analysis. 

The greatest success in this direction has been in two dimensions ^ . This is thanks to the following 
two facts. First, in two dimensions the scale invariance of a statistical system at its second order 
phase transition point, combined with rotational and translational invariance, implies full conformal 
invariance of the fixed point theory ||5|. Second, there are many completely integrable field theories 
which can be regarded as perturbations of these critical models off criticality Q]. Therefore the space of 
two dimensional field theories, particularly conformal and integrable ones, provide the best theoretical 
playground for addressing many of the issues mentioned above. 

Two dimensional conformal field theory, due to its central role in string theory and possible 
classification of all two dimensional critical phenomenal^, has been studied extensively in the past 

^For an example in two dimensions see[S|. 
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several yearsQ. Conformal field theories, besides describing critical statistical systems, give ultraviolet 
fixed points of the associated renormalizable quantum field theories. Though conformal invariance is 
enough to determine all correlation functions of the fixed point theories, these theories are by no means 
trivial. Unitary minimal models of BPZQ, the simplest of which correspond to the critical and the 
tri-critical Ising models, are all examples of such theories. Even though these theories admit free field 
representations, they are still non-trivial, in light of the fact that conformal symmetries only admit 
coordinate space realization^. 

Certain perturbations of minimal models are known to be completely integrable. The simplest 
example of this is the thermal perturbation of the two dimensional Ising model. A much less trivial 
example is the RG flow between tri-critical and critical Ising models. It is believed that a suitable 
integrable perturbation of the tri-critical Ising model has the critical Ising model as its IR fixed point |3|, 
^. Such flows, aside from being physically interesting in their own right, provide useful toy models 
for higher dimensional low energy effective field theories of completely solvable high energy theories. 
These provide ample reason for a need for a calculational scheme which naturally takes advantage of 
the existence of non-trivial yet solvable fixed point theories. Even though there has been some success 
in doing perturbative calculations near non-trivial fixed points of specific models up to low orders, they 
have not lead to a systematic calculational scheme up to now. 

Recently, expanding upon ideas originally put forth by Wilson [|l|, [l^, and expanded upon by 



Wegner 1 11 , 12| , a new framework for the analysis of the theory space of all renormalizable D dimensional 
euclidean quantum field theories has been introduced (]l3[-[p^]). Though this formalism is in principle 
non-perturbative, it provides a concrete program for perturbative analysis and calculation of the short- 
distance expansion of the correlation functions of a renormalizable field theory around a non-trivial 
fixed point. The first step in this program is to make a particular choice of coordinates for the infinite 
dimensional theory space such that all but a finite number of coordinates vanish for renormalizable 
theories. With this choice we can analyze the structure of the renormalization group equations by 
counting the scale dimensions of the coordinates [|^, 

The next step is to combine the renormalization group with the operator product expansion (OPE) 
and determine all non-analytic dependence on the renormalized parameters in the short-distance ex- 
pansion of the n-point correlation functions |l^]. At this point, all that remains for us to calculate are 
the RG invariant functions appearing in the solution of the RG equations 1 13]. Their evaluation requires 



an explicit calculation based on a concrete regularization scheme, which was provided by the introduc- 
tion of a variational formula for correlation functions in coordinate space [|^]. In this formalism the 
operator product expansion is introduced and incorporated in the definition of the variational formula 
in order to regulate all the ultraviolet divergences. It is not necessary to establish the validity of the 
OPE by extensive calculations [0] as it is already built into the formalism {i.e., OPE is assumed). 

The variational formula for the correlation functions realizes the derivative of correlation functions 
with respect to a parameter in terms of an insertion of a conjugate operator. The notion of conjugacy 
is the same as the one used in thermodynamics, in the sense that the expectation value of any such 
operator is defined to be the derivative of the free energy density of the system. The existence of such 
operators is an assumption. The variational formula is formulated in coordinate space, and involves 
divergent subtractions due to short-distance singularities of the renormalized theory. We also need to 
introduce finite counterterms to compensate for the arbitrariness of the subtraction scheme. The finite 
counterterms can be interpreted geometrically as connections in the theory space|16, 19|. 



See articles in ref.|^] also. 

For these models the standard Feynman techniques are not applicable. 
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The variational formula is very hard to derive from first principles and we can only demand its 
consistency with RG and the usual properties of derivatives. Requiring consistency with the RG equa- 
tions leads to a relationship between the anomalous dimensions of the composite operators and the 



unintegrable part of the operator product expansion [ 16 1. The curvature of the above connection ap- 
pears in this relationship naturally. The commutativity of mixed higher order derivatives (Maxwell 
relations) leads to an expression for the curvature of the connection in terms of a double integral over 
a finite domain. This provides us with a very natural setting for studying issues relating to general 
covariance and geometry in theory space[^, 19|. Upon applying the variational formula for correlation 



functions to the OPE we arrive at a variational formula for the operator product coefficients [|1 7]]. Re- 
peated application of these variational formulas allows for a perturbative calculation of the correlation 
functions. 

The large number of exact results in two dimensional conformal and integrable field theories mo- 
tivates the application of the above formalism to the study of the space of two dimensional quantum 
field theories. Since the above framework and its full conceptual and calculational power is not yet 
popular we decided to present a comprehensive discussion of the systematic calculational scheme which 
is contained in this formalism. In particular we illustrate the usefulness of this scheme in analyzing a 
historically significant model, the two dimensional Ising model. This model corresponds to the best 
understood non-trivial fixed point in the space of two dimensional field theories. It was solved exactly 
nearly half a century agoQ, and has served as an important toy model for the study of existence and 
definition of continuum field theories starting from lattice models Q. The correlation functions of the 
model are known exactly, thanks to its complete integrability at and away from the fixed point H, 
The short- and long-distance expansions of the correlation functions encode the operator content of 
the UV and IR fixed point theories. It is desirable to understand the structure of these expansions 
field-theoretically. Therefore, our interest in this model is twofold. First, since the presence of the spin 
field in the Ising model does not allow us to apply the traditional Feynman approach to perturbative 
calculations in field theory, the use of the new formalism is essential. Second, we can check the results of 
our calculations againt the available exact results. The crucial point is that this method is in principle 
equally powerful for many other interesting models for which exact solutions are not known. 

The purpose of the present paper is to discuss the formalism sketched above in a general renor- 
malizable field theory with a UV fixed point, and to extend the proof of infrared finiteness of the 
variational formula for the operator product coefficients to any order. This essentially amounts to 
the consistency of the variational formula with the analyticity of the OPE coefficients and leads to a 
manifestly finite variational formula for general operator product coefficients^ We therefore arrive at 
a systematic method for the perturbative calculation of OPE coefficients, and consequently correlation 
functions, in the vicinity of any non-trivial ultraviolet fixed point. As an application of our results 
we calculate the spin-spin and the spin-spin-energy density correlation functions in the massive Ising 
model perturbatively. 



This paper not only provides the conceptual and calculational details of our previous work[p6|. 



but also extends the calculational framework of ref. [|T^, 17| to higher orders. Our manifestly finite 



calculational scheme is a product of a complete examination of the infrared structure of the higher order 
variational formulas, and demanding consistency with locality principles. Calculation of non-trivial two- 
and three-point functions in the two dimensional Ising model within this scheme are explicit examples of 
the practicality of this formalism in addition to its conceptual clarity. These are the first field-theoretic 



*For the Onsager solution of the Ising model, see ref. and references therein. 

^These are the coefflcients appearing in the expansion of correlation functions in terms of expectation value of single 
operators. 
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calculations near a non-trivial fixed point ^and should, in principle, be generalizable to many other two 
dimensional models. The main obstacle for carrying out the same analysis for other models is finding 
a systematic way of calculating the necessary finite integrals, which are quite complicated in general. 
However, our success in finding a systematic technique for calculating the integrals for the Ising model 
correlation functions gives us hope that the same should be possible at least for other unitary minimal 
models. 



The organization of the paper is as follows. In section 2 we review the RG analysis of refs. |13]-[15| 



in general and in the two dimensional Ising model. In section 3 we review the variational formulas 
[]l6| , p^ ] and discuss their ultraviolet and infrared divergence structures in detail. We also introduce a 
variational formula for general operator product expansion coefficients in this section. The consistency 
of the variational formula with the differentiability of the OPE coefficients to a given order is shown 
to be equivalent to the IR finiteness of the variational formula to that order. Hence, the analyticity of 
the OPE coefficients implies that the variational formula is manifestly free of divergences to any order. 
In section 4 we analyze the UV structure of the variational formulas for the Ising model. Here we fix 
our choice of the parameters used in the remainder of the paper. In section 5 we apply the results 
of section 3 to the Ising model. We describe the structure of the n— order correction to the scaling 
behavior of the short-distance expansion of the spin-spin correlation function, and explain how the 
required integrals can be performed. For completeness, the details of the explicit calculation to third 
order [p6| are presented in appendix B. In section 6, as another non-trivial application of our results, 
we discuss the calculation of the spin-spin-energy density correlation function to first order. In section 
7 we present our conclusions. In appendix A some relevant facts about the operator content and the 
correlation functions of the critical Ising model are presented. 



2 Renormalization Group Analysis 
2.1 General Theory 

The Action of the Renormalization Group on the Theory Space 

In the Wilsonian approach to renormalization we consider the theory space (TS) of all regularized 
quantum field theories in D dimensional euclidean space Q. The theories in TS are regularized with an 
ultraviolet cutoff a which is assumed to be invariant under both rotations and translations. The cutoff 
may be taken to be a lattice unit as long as the continuum theory obtained in the limit of vanishing 
lattice spacing is rotationally and translationally invariant. Let Ti. denote a point in the theory space 
representing the action oi a D dimensional euclidean quantum field theory]]. 

The action of the renormalization group (RG), denoted by TZ{E), is parametrized by a real variable 
t and is defined as an integration over fluctuations between a and ae~^ followed by a dilation by a 
factor . Hence, the distance r, measured in units of the cutoff, is transformed to re~^, while the 
renormalization point r = 1 is fixed. Note that RG acts toward the infrared and differs from the 
standard definition by rescaling. We generally assume the existence of fixed points in TS with respect 
to the RG transformations. Let 7i* denote one such point defined with the property 

^Dotsenkoji^ has calculated the second order correction to the spin-spin correlation function, but, as he points out, 
the generalization of his analytical continuation method to other models seems formidable. 
^Equivalently 7i can be regarded as the hamiltonian of a D dimensional statistical system. 
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n{i)n* = n* (v^gr). (2.1) 



At the fixed point scale invariance implies 

(0n(Ari)---(^,„(Ar„))7^. = -^^(^— ^r:p^(0n(n) • • • (Ai„K)>w*^ (2.2) 

and we can assume that the operators (pi have non-negative scale dimensions xi, and are called relevant, 
marginal, or irrelevant if Xi are smaller than, equal to, or greater than D, respectively. 

The hamiltonian of a theory in the neighborhood of 7i* is obtained by adding a linear combination 
of local scalar^ operators to the fixed point hamiltonian. Let {(^j} denote a complete basis of local 
scalar operators which are infinite in number. The generalized mass parameters {^*}, introduced by 
Wilson, are the coefficients of these local scalar operators and provide local coordinates for the TS. 
Hence, TS is infinite dimensional. 

Generally, starting from any hamiltonian Ti in the neighborhood of H*, RG transformations 7?.(^) 
generate a trajectory of hamiltonians TZ{t}7i parametrized by The range of £ can be extended in both 
directions as long as the trajectory stays within the theory space. Not all trajectories can be extended 
without limit. A physically interesting subspace of TS is the space of renormalized trajectories, denoted 
by 5(00), consisting of all trajectories which can be traced back to ^ = —00 without going out of the 
theory space. ^ Renormalizable theories are typically characterized by a finite number of renormalized 
parameters and we should therefore be able to parametrize S'(oo) by a finite number of parameters 
forming a finite dimensional subspace of the theory space, accessible asymptotically in the £ — > 00 limit 
of the RG flows. Let us explain these points more carefully ||l|, [2l]] . 

Consider the lattice models parametrized by couplings |^ where i G T with T an (in general infi- 
nite) index set. Assume that the model undergoes a second order phase transition at k'^^-'* parametrizing 
the critical hamiltonian Ti-iJ^^^^^. At the critical point all parameters have well-defined scale dimensions. 
Hi = D — Xi, and we assume that only a finite number of them have positive scale dimensions. Starting 
from an almost critical theory, under RG transformations the effect of the couplings with negative scale 
dimensions gets washed out. Therefore in the infrared limit the action of the renormalization group 
singles out a finite dimensional subspace of the theory space corresponding to a renormalized trajec- 
tory, i.e., universality. The dimension of this subspace is equal to the number of relevant operators 
conjugate to couplings with positive scale dimensions. 

When the continuum theory exists, with the proper choice of coordinates (which will be described 
shortly) W^^l^j^g flows to the origin of this subspace corresponding to H* which is the UV fixed point of 
a renormalized field theory originating from it. In practice after reaching the renormalized trajectory 
in the infrared limit we can trace back to 7i* along the renormalized trajectory]^. Therefore these 
theories are defined up to infinitesimal distances and on the trajectory the parameter a plays the role 
of a renormalization point. In the process of defining the continuum limit in this way we are led 
to considering a finite number of renormalized parameters, g^, characterizing this trajectory, which 
are analytic functions of the couplings k*. Therefore the goal would be to find a particular choice of 

^Thanks to the assumption of rotational invariance we need to only consider scalar operators. 

^This is equivalent to the assumption of the existence of a UV fixed point at the origin of renormalized trajectories. 
^"This is one particular choice for Wilson's generalized mass parameters. 

^^Since proving the existence of UV fixed points is in general very difficult, we assume it. In principle these fixed point 
theories are fully interactive. 
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renormalized parameters, which allows for the study of the non- analytic dependence on the generalized 
mass parameters in the short-distance expansion of the correlation functions. In the next subsection 
we will see how this happens for the Ising model. 

While Wilson's generalized mass parameters provide local coordinates for the infinite dimensional 
theory space, they are not necessarily the most suitable choice for a smooth parametrization of the 
renormalized trajectories. First, since we only need a finite number of renormalized parameters to 
parameterize the renormalized trajectories, we would expect that there must exist a new set of coordi- 
nates for TS such that all but a finite number of them vanish on ^(oo). Second, it would be desirable 
to be able to analyze the structure of the RG equations for these new coordinates by counting scale 
dimensions only. The advantage of these new coordinates would be that their scale dimensions are 
additively conserved under RG flows. This, in addition to the analyticity requirement of the RG equa- 
tions implied by the locality of the theory, restricts the form of the RG equations to contain mixing 
with parameters of lower or equal scale dimensions only. 

Based on Wegner's scaling fields pi , it has been shown that the analytic transformation of interest 



to the new coordinates g^{i = 1 • • • A^) exists |13]. Under this change of coordinates, which is assumed 
to be analytic, our basis of local operators transforms like a vector. Denote the new basis of composite 
operators by {(j)i}g- In terms of these new coordinates, the operators (pi have well-defined short-distance 
singularities only for the theories on S{oo). This allows for a systematic determination of the corrections 
to the scaling behavior of the short-distance expansion of the correlation functions. Note that if we 
consider general theories not restricted to ^(oo), the correlation functions contain arbitrarily high 
singularities, even in terms of the new coordinates. 

The coordinates of the theory space can be chosen such that the action of the renormalization group 
on the renormalized trajectories is 



■^/ = P\9)^y^9'+ E ^/^iV, (2.3) 



where / stands for ii ■■■ in, yi = yh + ■ ■ ■ + yi„ and = g^^ ■ ■ ■ g^" ■ The coefficients Pij are independent 
of g. Note that g = is kept invariant under (2.3). In ref. [13| it was shown that with this choice of 



coordinates the short distance limit reduces to the limit of ^f* ^ 0. 



Vector Bundle of Composite Operators 

Each point g in the space of renormalizable theories is a renormalized field theory which has an 
infinite number of linearly independent composite fields. {4>a}g denotes a complete basis of composite 
fields at g. The space of composite fields is a linear space, and therefore forms an infinite dimensional 
vector bundle over 5(oo)|jl^. The choice of the local basis is by no means unique, and two bases {4>a}g 
and {4>'a}g are related by an invertible infinite dimensional matrix N{g) via 



0'a = (A^(5))aVfe. (2.4) 

The physical quantities, the correlation functions of n composite fields {(paiifi) ' ■ ■ 4>a„i'>'n))g, form 
rank-n tensors on the theory space. 

Among the infinite number of composite fields, composite scalar fields stand out; they are the 
fields Oi conjugate to the parameters g^. The notion of conjugacy should be familiar from thermody- 
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namics. The derivative of the free energy density F{g) with respect to is given by the expectation 
value of 



= ^Fig). (2.5) 
The RG action (|2.3| ) combined with canonical scaling of the free energy density 

^F(<7) = DF{g), (2.6) 
and (|2.5|) induces an RG equation for composite operators O,: 

d 83^ 

_0. = DO. - I-O,. (2,7) 

The particular choice of parameters {g*} restricts the form of RG equations for the general composite 
fields; scale dimensions are additively conserved, and under RG any field can only mix with operators 
of equal or smaller scale dimensions, with mixing coefficients analytic in the parameters. Therefore, 
introducing an infinite column vector $ of composite operators, we can write 

^<^> = r(5)'J>, (2.8) 

where the matrix 

^a\9) = ^aSa' + 0{g) (2.9) 

gives the scale dimensions Xa and the anomalous dimensions of renormalized fields. The precise meaning 
of the above RG equations is given by the following RG equations for the correlation functions: 



E ^a,\9){ct>a, in) ■ ■ ■ Urk) ■ ■ ■ <^a,. {rn))g. (2-10) 



k=l 



Operator Product Expansion 

The operator product expansion (OPE) plays a central role in the analysis of the short-distance 
expansion of correlation functions p!o|]. Wegner combined OPE with the RG equations to ascribe all 
non-analytic dependence on the scaling fields to the vacuum expectation values of single composite 
operators |[T^. When r defined by = X]j<j(^j ~ is small compared with the correlation length, 
we can write 



^^In order to make the notion of conjugacy in field theories precise, we also need to specify how the derivative of a 
correlation function is realized. The variational formula (VF) in section 3 is one such realization. 
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'/'ai (n) • • • (t>a„ {rn) = Ca^-a„^ {n, ■ ■ ■ , r^, g)h{^), 



(2.11) 



where Ca^-.-aJ' is invariant under translation. The vahdity of this equation as an operator equation 
requires the sum over b to run over ah local operators with or without vanishing spins. This is 
understood as an asymptotic expansion of the correlation functions 



Kin) ■ ■ ■ 0a„(r„)0d,(ei) • • • (l^dM))9 = Ca,-aJ'{ri ■ ■ ■ r^, g) (0b(O)0d,(ei) • • • (t^d,{ik))9. (2.12) 



The fundamental assumption for the operator product expansion is that the coefficients Cai...aJ'{g) 
depend on the parameters only analytically. This is equivalent to the local nature of the operator 
product expansion: short-distance physics alone guarantees ( 2.121) to be valid. 



Using the RG equation (2^) we can derive the RG equation for the coefficient functions 



,/,Cai---a„ (TIi ' ' ' I'^n] g) — Cai - a„ ^d (ff) ~l~ '^aj {g)Cai---ai-idai+i --a 



(2.13) 



i=l 



where the matrix T{g) is given by (2.9). 



Solution of RG Equations 



In order to solve RG equations (2.10) and ( 2.13| ) we define the running parameters g{t,g) as the 
solutions of RG equations 



d_ 

de 



g\e-g) = f3'igii;9)), 



(2.14) 



with the initial conditions g^{0;g) = g^. We will suppress these initial conditions from now on, and we 
will denote g^{i',g) simply by g^{f). We can then construct RG invariant quantities, ^*(lnr). The 
RG invariance follows from the fact that the dependence on the change of the coordinate r under the 
RG is cancelled by the dependence on the change of initial parameters g under the RG. 

We also need to introduce the matrix G{r;g) that satisfies 



d_ 



Gir;g)=Tig)Gir;g), 



(2.15) 



with the initial condition G{l;g) = 1. The matrix G{r;g) is easily calculable for any field theory with 



at most one marginal parameter. Using this matrix the solution to equation ( 2.13 ) is 



Ca,...aJ{nr ■ ■ ,rn;g) = Ga,'Hr;g) ■ ■ ■ GajHr;g)Ht,...K'{^^^^rg{^^^^^ (2.16) 



where 
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is an analytic function of g, due to the analyticity assumptions for Cai---a„'^{g)- Therefore the solution 



of (|2l0D is 



{cl>aAn)---^aArn))g = Ga,'Hr-,9)---Ga^^{r-g)H,,...,/[^^^ (2.18) 

and all the non-analytic dependence on g's in the correlation functions comes from the expectation 
values of single operators through the non-analyticity of the free energy density. 

Since the short-distance limit on renormalized trajectories is equivalent to the limit ^ 0, we 
obtain the short distance behavior 

{<Pa, in) (r„)), ^ {cPa, (n) • • • 0a„ (r„))w* . (2.19) 

which is consistent with our previous comment that the operators (pi have well-defined short distance 
singularities only for the theories on S{oo). 

2.2 The Two Dimensional Ising Model 

The Ising model is defined on a periodic lattice by the hamiltonian 

N 

n+l}, (2.20) 

m,n=0 

where am,n = ±1. This theory undergoes a second order phase transition at K = Kc given by 

sinh {2Kc) = 1. (2.21) 

The modern renormalization group language offers a qualitative description of continuous phase 
transitions. Let us describe the phase transition in the Ising model in this language. TLiatUce defines 
a one parameter family of theories in the theory space. The action of the renormalization group, in 
the infrared limit, singles out a renormalized trajectory in the theory space such that theories with 
K < Kc and K > flow to two different tails of this line. The critical theory 'HiJ^nce^ ^'^^ ^ ~ 
flows to Tl* (the origin of the coordinates on the renormalized trajectories) along a critical line. On 
the renormalized trajectory, the vicinity of the origin characterizes the short-distance behavior of the 
scaling limit of the lattice theory. This limit is defined by letting e = 4:{Kc — K) — > 0, but keeping 
the product t = eR {R is the distance in lattice units) finite. The renormalized mass m should be 
proportional to e up to first order. The complete dependence of m on e is fixed by demanding a simple 
RG equation for m. 
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In a statistical mechanical system near its second order phase transition point, where the correlation 
length diverges, information about the details of the lattice is washed out. If the above scaling limit 
exists we have a well-defined quantum field theory all the way up to zero distance. This fact is 
usually stated as the existence of an unstable UV fixed point which can be reached backward along 
the renormalized trajectory, using the RG equations for the renormalized parameters. We recover 
rotational invariance in the continuum limit, and the theory is equivalent to a theory of a free massive 
Majorana fermion. The theory is parametrized by the cosmological constant gi and the mass m, 
which are conjugate to the operators 1 and Om (corresponding to the energy density operator iipijj): 
respectively. These parameters satisfy the following RG equations 

^91 = 2gi + ^(3im^ (2.22) 
= m. (2.23) 



There are two RG invariant quantities 



In this case (|2.7D reduces to 



gi (In r ) = r'^gi + ^/3im^r^ In r (2.24) 
m.(\nr) = mr. (2.25) 



4Cm = 0„-/?iml. (2.26) 



Taking the expectation value of both sides and using ( 2.23| ) we have 



{0^)m = -l3imln{mc), (2.27) 
where c is a constant of integration. Equations ( ^.27 ) and (|2.5| ) give 



F{gi,m) = 51 - i/5im2 In(mc) + ^Pim\ (2.28) 



This was one of the main results of ref.[14 



Since gi only contributes an additive constant to the free energy density, none of the correlation 
functions of local operators depend on it. In this case, we can effectively regard ^(cxd) as a line 
parametrized by a real parameter m. This is the line corresponding to the theory of a free massive 
Majorana fermion mentioned above. Therefore from now on all of our physical quantities will depend 
on m only. 

The RG analysis of one point functions of more general composite operators in the Ising model is 
quite straightforward. The first reason for this, as shown in appendix A, is that the couplings to all 
the members of class a can be turned off, due to the locality of the RG transformations and the special 



OPE structure of the operators of the critical Ising model (see A^). The second reason is that, due 
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to the additive conservation of the scale dimension of operators under RG transformations, the matrix 
r(m) has the following structure 



Ta^m) = Xa5a^ + ja^m"" iff X, - G Z+ (2.29a) 

= XaSa'^ otherwise. (2.29b) 



Note that ( p. 29a ) is particularly simple since there is no marginal {i.e., dimension 2) scalar operator 



in this theory. Using ( |2.23| ) and (2.29a) we can write (2.8) as 



^[m-^(<l>)™] = 7[m-^(<I>)„^], (2.30) 

where we have introduced the infinite dimensional diagonal matrix X = diag[xi]. 

Now since we are considering a unitary model, the lowest dimension operator is the identity oper- 
ator with scale dimension zero. Arrange the column vector $ such that operators with higher scale 
dimensions are lower in the column. With this choice the matrix 7 is lower triangular. Therefore upon 
repeated differentiation of the above we get a non-zero result 

^[m-^'^i^aU = la'^lt^ ■ ■■lKJ"[m-'''H<^bJm], (2.31) 

if and only if a > 61 > • • • > 6„. We therefore see that the case 4>b„ = 1 corresponds to the largest 
value n can have. Using equation ( 2.23| ) we see that n is the highest power of In(mc) and is equal to 



the number of scale dimensions smaller than Xa, where at least one non-derivative scalar composite 
operator exists. We therefore have 



^7a'^76,'^---7fe„_i'"[ln(mc)]" + ---. (2.32) 



Recall that we can change our basis of local operators, which is equivalent to an analytic change 
of coordinates in the neighborhood of renormalized trajectories. In other words, this transformation 
should not spoil the mass independence of the RG equations 

k = (t>a+ E nJ'nf^-'"'ct>i, = Na\m)(t>i„ (2.33) 

b:Xa>Xi, 

such that with 7 = N{m)'yN~^{m) we can write: 

{<Pa)m = m^'^^^ja'^^"^ ■ • • 7fe„_, [ln(mc)]^ (2.34) 



n 



This transformation exists if a reduced matrix, obtained from 7 upon eliminating the entries due to 1, 
is diagonalizable. 

For the RG analysis of more complicated correlation functions we need to know the matrix G{r; m). 
In fact we can solve ( 2.15| ) very easily to get 
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G{r; m) 



X 

m r 



[mr 



-X 



(2.35) 



where 7a'' are non-vanishing only if Xa — x;, is a positive integer. With the exphcit form of G{r;m) at 
hand, the correlation functions ( ^.18 ) are determined up to RG invariant functions. The remainder of 
this paper is concerned with perturbative calculation of these RG invariant functions. 



3 Perturbation Theory for OPE Coefficients 



As we have seen in equation ( p. 16 ), after separating out the non-analyticities of correlation functions 



encoded in expectation values of single operators, we remain with the OPE coefficients which are 
analytic functions of parameters. Renormalization group analysis fixes these analytic functions up 
to RG invariant combinations (2.17). Determination of these functions requires explicit perturbative 



calculations, based on a concrete regularization scheme. Note that, since the OPE coefficients are 
universal, they do not depend on how we calculate them. 

On the other hand, conjugacy in thermodynamics and statistical mechanics is an exact notion, 
irrespective of any approximation technique used in practical calculations. The thermodynamical def- 
inition of conjugacy is that the derivative of a correlation function with respect to a parameter is 
given by an insertion of the volume integral of the corresponding conjugate operator. As we already 
discussed, well-defined renormalizable field theories are defined as the scaling limit of statistical sys- 
tems. In taking this limit, physical short-distance singularities appear, which prevent us from naively 
borrowing the definition of conjugacy as given in thermodynamics. In order to extend this definition 
to field theories, the notion of singular operator product expansion is introduced to cure the physical 
ultraviolet singularities mentioned above. Moreover, the full operator product expansion (both singular 
and integrable) need to be introduced to incorporate locality into this formalism. 

The variational formula for correlation functions gives a realization of the derivative with respect to 
a renormalized parameter which satisfies the requirements imposed by the field-theoretical conjugacy 



relation|16]. Consistency of this variational formula with the locality of the operator product expansion 
amounts to the infrared finiteness of the variational formula for the operator coefficients. The latter 
variational formula has been introduced and proven to be infrared finite in ref.[0]. We provide the 
details of the proof of the infrared finiteness to all orders. As a by-product we arrive at a manifestly 
finite calculational scheme for the perturbative calculation of the OPE coefficients. In the following 
sections we discuss these issues in detail. 



3.1 Variational Formula for Correlation Functions 



In order to characterize the short-distance singularities of the renormalized correlation functions, 
whose volume integral will be needed shortly, and to choose the suitable subtractions to make these 
integrals UV finite, consider the following singular OPE coefficients 
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0,{T)U^) = Ctl\r-g)U^)+0 



(3.1) 



Here, we assume the average over the orientation of the coordinate vector has been taken, such 
that the right hand side includes only the scalar composite operators. Also note that only composite 
fields (ph whose scale dimension are at most Xi + Xa — D should be kept; we only need to consider the 
unintegrable part of the OPE coefficients for calculating the subtractions. 

The divergences arising from the singular OPE coefficients are physical short-distance singularities 
and they exist even after renormalization. By using these singular OPE coefficients, the following 
variational formula has been introduced|l(;] 



-5j(0ai(n) • • • (Panirn))g = hm 



d^r' {Ot{r')<t>aAri)---<t>aSrn)), 



|r'-rfc|>e 



+ ^iaj'ig-. e)(0ai (n) • • • 0fe(rfc) • • • (t)aSrn))c 



k=l 



(3.2) 



where di = the definition of connected operators is 



Of(r)^0,(r)-(0,)„ 
and the definition of the subtractions, in the matrix notation, is 



(3.3) 



l>r>e 



d''rC['\r-g)+Ci{g). 



Here 1 is the renormalization point, fixed under the renormalization group. 



(3.4) 



The arbitrariness involved in the subtraction of the divergent counterterms is compensated for 
by the introduction of finite counterterms Ci{g). These finite counterterms, in addition to the above 
mentioned physically important role, have a natural geometrical meaning in the theory space. The 
correlation functions are tensor fields on the space of renormalized trajectories, and therefore their 
partial derivatives are not covariant quantities. To be able to define a covariant derivative we need to 
have a connection on the space of renormalized trajectories in theory space. In other words, covariance 
of the VF demands the introduction of finite counterterms Ci{g) which transform as a connection 



c,{g) ^ c[{g) = N{g){d^ + c^{g))N-\g), (3.5) 

under the change of the basis given by equation ( |2.4| ) . This is how connection q {g) arises naturally in 
the space of renormalized trajectories[^]. 

Consistency Conditions of VF with RG Equations and Maxwell Relations 

It is crucial to notice that the above VF is assumed to be valid for finite values of the parameters 5', 
i.e., to be valid beyond perturbation theory. This is a reflection of the non-perturbative nature of the 
conjugacy relation even in field theory. Though the definition of the variational formula is inspired by 
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its statistical mechanical counterpart, the relation between a lattice model and the corresponding field 
theory is non-trivial. Consequently, it is difficult to derive the VF from first principles. Instead, if the 
existence of the conjugate fields Oi, satisfying the conjugacy relation (|2.5D and VF (3.2) is assumed. 



then the consistency of VF can be checked. By demanding consistency with the RG equations (^.3|), 
(pJl) , and (pj) we arrive at the following remarkable relation |]lf 



voliS^-') Ht\g) = -^ng) + Hg),^ig)]+P^ig)n,,{g), (3.6) 



where 



^ig)^r{g)+(3\g)c,{g), (3.7) 



is a covariant tensor field, 



Hl^\g)^Ct'(r = l;g), (3. 



and the curvature Qij is defined by 



^ijig) = diCj{g) - djCi{g) + [ci{g), Cj{g)]. (3.9) 

The curvature could be determined by demanding the consistency of the VF, when used recursively 
to evaluate higher order derivatives. This means that we want the VF to satisfy the Maxwell relations: 

didj = (3.10) 

Imposing these conditions on the second order VF for correlation functions gives an expression for the 
curvature in terms of a double integral over a finite domain[16|: 



a,(5)(^)9 = / d^'r F.P. / d^r' (Oj(r) {Ofir') - Q(r')) $(0) - (i - j))„ (3.11) 

Jl>r Jl>r' 

where F.P. denotes taking the integrable part with respect to r. The commutativity condition for more 
than two derivatives reduces to the Maxwell relations ( p. 10 ) and does not give any new conditions. 



Consistency condition (|3.6| ) gives a geometrical expression for the singular part of the OPE coef- 
ficients; the vector field P^{g) (beta functions), the rank- two tensor field ^{g) (scale and anomalous 
dimensions), and the connection Ci{g) (finite counterterms) , determine the unintegrable part of the 
short-distance singularities of the theory. These are analogues of the well-known relation between 
the renormalization constants, the beta functions, and the anomalous dimensions in the dimensional 



regularization with minimal subtraction scheme] 22]. 
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3.2 Variational Formula for OPE CoefRcients 



Even though the above VF for correlation functions is assumed to be true non-perturbatively 
(for arbitrary values of g), it is hard to use in practical calculations. However it is naturally suited 
for calculating off-critical correlation functions around a non-trivial fixed point by using only the 
information about the correlation functions at the fixed point {g = 0), where they are easier to calculate 
because of invariance under scale and sometimes larger symmetry group transformations. This is 
nothing other than perturbation theory around a UV fixed point, and therefore requires the introduction 
of a realization of derivatives with respect to the parameters of the theory, valid particularly at the 
point g = 00. 

The non-analytic dependence of the correlation functions on the parameters shows up in the VF 
as IR divergences at g = 0. But by virtue of ( p. 18 ) all the non- analytic dependence of correlation 
functions are encoded in one-point functions. Therefore we have to look for a VF for OPE coefficients 
which, by assumption, are analytic in parameters and perturbatively calculable. This way of doing 
perturbative calculations is in the spirit of ref. [10| and more general than the standard diagrammatic 
way. 



In ref. [17|, using the VF of correlation functions, a VF formula for the OPE coefficients in 
the product of two operators was given. Along the same line using the VF (^]^) for the derivatives 
of correlation functions on both side of OPE ( 2.11| ), we obtain the following VF for general OPE 
coefficient^ 



-9iCai...a„ {ri, - ■ ■ ,rn;g){(l)b)g = 



lim 



n 



k=l 



(3.12) 



To be consistent with the differentiability of OPE coefficients, the above VF should be free of IR 
divergences. This consistency is automatic because from OPE (|2.1l]) , for any R' > maxjri, • • • r^}, we 
have 



d^r' (Of (/) Ua, (n) • • • (r„) - Ca,...aJ{n, • • • , r„)</),(0)) ), 

r'>R' ^ ' 



(3.13) 



So if we consider B! as the IR cutoff for equation ( 3.12| ), the integral is independent of the cutoff. On 
the other hand, since the integral is performed over a finite domain, it is free of infrared divergences as 
well. As expected physically, this guarantees that UV physics completely determines the structure of 
the OPE coefficients; B! can be chosen as small as, but not equal to the largest argument of the OPE 
coefficients. 



Usually the correlation functions are known not only for the fixed point but also for some special path out of the fixed 
point. In these cases we can use VF on the trivial path to calculate the correlation functions in its vicinity perturbatively 
(e.j., In the massive ((;/)*)4 theory, since all correlation functions are known for the free massive theory, the VF for A = 
could be used to peturbatively calculate the OPE coefficients in the vicinity of the free massive path |17[). 

^*For the rest of this section all formulas are for a finite but for simplicity we suppress some of the ^-dependences. 
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The VF ( p.l2D easily generalizes to higher-point functions. We have the following first order VF 
for general OPE coefficients when max{ri ■ • ■ r„} < ii' < min{si • • ■ s;} 



diCa^...aJ{ri, ■ ■ • ,r„; £/)((/){, (0)(/)di(si) 

i?'>r'>E 



lim 

e-*0 



k=l 



■(t>di{si))g = 

■ 4>aArn) - Ca^...aJ{ri, ■ ■ • ,r„)4(0)) 4>dAsi) ■ ■ ■ 4>di{si))g 
, rn)-Ca,...aJ{ri, • • • , rn)Aj{e)') {MO)<Pdi (si) ---(pdi {Sl))g 



(3.14) 



which is manifestly free of any divergences (both UV and IR). 

Higher Order VF and Consistency with Analyticity of OPE Coefficients 

We can use the VF recursively to evaluate higher order derivatives of OPE coefficients. The consis- 
tency of the higher order VF with the commutativity of the partial derivatives {i.e., Maxwell relations) 
has already been checked and resulted in an expression for the curvature. This type of consistency is 
related to the UV behavior of the variational formula; the curvature is in terms of double integral over 
a finite domain. But we have another consistency check: the VF should be consistent with the ana- 
lyticity assumption of the OPE coefficients. This is related to the IR behavior of VF. In the following 
we will show that, like the first order VF, this consistency can be made automatic for higher orders 
recursively. 

To avoid too many indices and to see the structure of higher order VF in a more transparent way, 
let us consider the OPE coefficients for two fields in the matrix form. Later we will only use this higher 
order VF for the Ising model. Consider the full OPE for two composite fields 



Consider the following two examples of first order, manifestly finite VFs for Ca- 



(3.15) 



-diCa{r;g){<^)g = Vim 



R'>r'>€ 



+ (^[Ai{e),Ca{r)]+Aj{e)Cdir)) (<!>) 



(3.16) 



-5,C7a(r;5)(<I>(0)</.rf(s))g = lim 



d^r' {Of{r') (Mr) - Ca{r)) mUs))g 

L " R'>r'>e 



+ {[Me),Ca{r)]+Ai/{e)Cd{r)) {^{0)Us)), 



(3.17) 



where r < R' < s. 
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The second order VF could be derived by taking the derivative of VF ( p.l6|) and using the VF for 
correlation functions (^.2|). At the end we have 



djd,Ca{r;g){^)g = lim 



d r hm 

l'--'-|>^ n^O 

- R'>r'>e 



J?">r">T) 



+ 



(d,Ca[r) + \A,{ri).Ca{r)\^A,a\ri)C,{T)) (Of (r')a>(0)), + V(^)('^^(O(0a(r) - ^(r)) $(0)), 
+^,,''(r?)(Of(r')(<^.(r) - C,(r))$(0))3 + A,{r^){pt{r')[Ur) " C,(r)) <I>(0)),) 



X lim 

?7^0 



R">r">ri 



(3.18) 



where r < R' < R" . It is important to note that in the above formula, and all higher order VFs, 
whenever two or more connected operators are inserted, partial connectedness with respect to all of 
the inserted operators is implied. 

The R' independence of VF ( |3.18| ) is guaranteed since our starting point, the first order VF ( 3.16| ), 
is already shown to be independent of R' as a direct consequence of the OPE itself {i.e., zero-th order 
VF). In the same spirit we can use the first order VF contracted with two point functions ( 3.17| ), to 
find the following identity 



lim 



7' — r I > e 
L " R'>r'>€ 



d r 



r">R" 



+ 



d^r" {0'^^{r")Ot{r')(<t>a{r) " C,(r)) $(0))^ 
(d,Ca{r) + [A{e),Ca{r)]+ Aj{e)C,{r)) f d^'r" (O,^(r")$(0)), 

^ ^ Jr">R" 



0. (3.19) 



This implies the R" independence of the VF ( 3.18| ). So we end up with a second order VF which is 
independent of IR cutoffs and its integrals are over finite ranges, and therefore manifestly finite. 

The above arguments for IR finiteness can be generalized, in a straightforward manner, to higher 
order VFs for general OPE coefficients recursively. Generally for n— order VF of OPE coefficients we 
end up with n nested volume integrals with increasing cutoffs from the outer integrals to the inner ones. 
The resulting expression is independent of the cutoffs, the smallest of which may be chosen as small 
as, but not equal to the largest arguments of the OPE coefficient. This proves the differentiability, to 
any order, of all OPE coefficients starting from the VF. In another words, this shows the consistency 
of the VF with the analyticity of OPE coefficients with imposing a mild but calculationally crucial 
new condition on the nestedness of the IR cutoffs; this variational formula respects the locality of the 
theory. 



3.3 Systematic Perturbation Theory around a Fixed Point 



Higher order VFs such as ( |3.18| ) give us a compact equation for the derivatives of the OPE coeffi- 
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cients. But, we need more refinements to develop a systematic method for perturbative calculations^]. 
First, we must show how to truncate the infinite sums in Second, since calculating the integrals 

with finite IR cutoffs is generally difficult, we will remove the cutoffs carefully. In this section, the 
above two refinements will be made simultaneously. It then remains to show that the VFs provide us 
with sufficient conditions for calculating the relevant quantities. We will explicitly show this to be true 
for the IM in section 5. 



Truncation of VF for OPE coefficients at the Fixed Point 

On physical grounds, we expect that, in the perturbative calculations around the point g = 0, the 
contribution of operators with high enough scale dimensions to variational formula for OPE coefficients 
should die off in the g limit. To show that VFs for OPE coefficients, like equation ( 3.16| ), have 
this property we need to go back to their derivations. To keep the presentation simple, let us again 
consider the case of the second derivative of OPE for two fields. First note that scale covariance 
at the fixed point gives 



{(l)a) =0 if <^a ^ 1. (3.20) 

Here and in what follows we use the convention that fixed point correlation functions are denoted 
without any subscripts. Taking derivatives of ( |3.15| ) we have 



^J^iCab''{r;g)mg = ^j^i{Mr)M0))9-^^Cab''{r;g)^j{4>d)g-^,Cab''{r;^ 

(3.21) 

Because of the basic analyticity assumption for the OPE coefficients the left hand side of equation 
(|3.21| ) is analytic in g except for {4>d)g- So, using ( |3.20| ) the right hand side is regular at 5 = 0; 
all divergences of the terms on the right hand side exactly cancel one another. Hence, ( |3.21| ) is well 
defined even at the fixed point. If we use the VF for the derivatives in equations ( 3.21| ), we will get the 



second order VF ( ^.18 ) with the infinite sums over operators with arbitrarily high scale dimensions. 



Instead, let us see what we can say about the derivatives of one point functions before using the VF. 
Translational invariance and dimensional analysis considerations imply 



d^{<Pa) = if xa>yi, (3.22) 

djdi{(l)a) = if Xa > yi + yj. (3.23) 



Using the VF for correlation functions, equations ( 3.22 ) and ( |3.23| ) give us a hierarchy of identities 



involving single and double volume integrals over correlation functions. These integrals do not have 
any IR cutoffs. Each of these integrals could be split into a main part with an IR cutoff and a correction 
part which is over the region beyond the IR cutoff. We can use these identities to truncate the VF 
( 3.18| ) such that only a finite number of terms remain. These identities also provide some restrictions 



on the choice of finite counterterms [Pq] . 

If we consider well-behaved field theories in which we only have a finite number of composite 
operators with scale dimensions less than any given finite number, then we are only left with a finite 

^^The notion of conjugacy in theories with conjugate operator(s) which are singular at y = is much more intricate |16[. 
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number of terms in the sums and all other terms could be written in terms of these corrections. Because 
of the dimensional restrictions and the connectedness of the insertions, the corrections vanish at g = 
when we take the nested infinite limit of the cutoffs (from the inner to the outer ones). After all these 
steps we end up with the following VF at the fixed point 



dACab\r)= lim hm \Ay'+y.--^-»)MijJ-,—) 



SVi+yj-^dn) 



+ r 



^'^^'>d,C,,'^ir)SlX(^) - {{i,R') - [, (3.24) 



where the sums on the repeated indices are restricted by 



Xd,, < Vi + Vj, 

and the main term and the two types of IR subtraction integrals are defined as follows: 



(3.25) 
(3.26) 



Mij,abiU', U") 



C-0 



U">U'l>i 



U'>u'>C 

-dj[A.,a%rOC,hHr)+Aih%rOCae\r)) + {A,a%rC)Cj^ (r) + A.^rOCae"^ (r)) 

>, (3.27) 



X hm / d^u" {O^{ru")cl>d,m+MHr0 

\Ju">u">i 



Si^{U') = r(^»+^'') lim^ 



l«'-l|>C 
- " (7'>u'>C 



(3.28) 



Sj^AU', U") = r(^»+^i+^-+^i') hm 



U'>u'>C 



^,_^^,^d^u' Ihn ( d^u" {O'^{ru"mru')Mr)M0)) 



u'l-u'\>i 

u">u">e, 



+Aji'{rOCab''{r)mru')M0)) + Cab''(r)Ajd%rC){Of{ru')M0)) 



-Cab\r)djAi/{rC) - Cab''{r)Aif{rC) lim / d^'u" (Oj(rtx")0e(O)) + A,e\C) 

^-♦0 \Ju">n">l; 



>, (3.29) 



where the range for the sums on the upper indices of the subtractions (e in the above equations) is 
restricted by the dimensional constraint in the singular OPEs (p.l|). Therefore, (3.24) is a finite (both 
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UV and IR) and truncated variational formula at the fixed point. After taking the limits as specified 
above, we do not have any cutoffs left, neither UV nor IR. The generalization of the above construction 
to higher orders is straightforward. It is, however, important to notice that the order of the limits is 
crucial; they do not commute. In addition, it is important that we only need to calculate the integrals 
for infinitely large IR and small UV cutoffs. The case of the Ising model is an instructive example of 
these points. 

Calculational Scheme 

The generality of the finiteness proof, presented above, assures the validity of the following shortcut 
for writing a truncated variational formula at the fixed point. We truncate the formal expression 
obtained by repeated differentiation of the OPE and its evaluation at the fixed point. Without loss of 
generality, consider the case of second order perturbation. From equations ( p. 21 ), ( 3.22| ), and ( 3.23| ) 
we find 



djdiCab\r) = djdi{Mr)(f>bm - djCj'{r)di{4>d^) - d,Cj' {r)dj{4>d,) - CjHr)dA{4>d,,) ■ (3.30) 

Now, we replace the derivatives of the correlation functions with their respective variational formulas in 
the following special way. Every di represents an integration over the position of the inserted operator 
Oj. The range of integration should be bounded by the same IR cutoff, respecting the nestedness of 



the cutoffs in the main term. So we arrive at the second order finite and truncated VF ( p. 241) , directly. 

Then, after having the VFs up to the desired order, we must go through the following steps. First, 
we should calculate the UV subtractions ( |3.4[) in terms of beta functions, anomalous dimensions, and 
finite counterterms using the consistency relation ( p.6| ) and the solution of the RG equation ( 2.16| ) 



for the singular OPE coefficients! 17]. Second, we must calculate all the necessary critical correlation 
functions. This can be hard for non-trivial fixed points, and is usually only possible for fixed point 
theories invariant under a large symmetry group. Third, we can calculate the volume integrals over 
these correlation functions with infinitely large IR and small UV cutoffs. This is usually the hardest 
step. We explicitly implement these steps for the spin-spin correlation function in the Ising model to all 
orders of perturbation theory. Finally, substituting the above relations into the derived VFs, we end up 
with a set of algebraic equations in terms of the scaling OPE coefficients, their corrections, anomalous 
dimensions of the composite operators, and the finite counterterms. After choosing a convention for 
off-critical composite operators by fixing the finite counterterms, up to the restrictions given by the 



explicit expression for the curvature ( ^.11 ), we expect to have sufficient information to solve these 



equations. A perturbative expansion of the short-distance expansion of correlation functions follows. 
This step has been carried out up to third order in our previous work ||2^ . 



Let us close this section by addressing a subtle point in using the VF at a fixed point. There are 
cases where one or more derivatives of some one point functions are non-analytic at the fixed point, 
but the naive 5 — > limit of the terms corresponding to these derivatives inside the VF give zero due 
to the critical OPE algebra. The reason for this apparent inconsistency is that, in order to capture the 
correct non-analyticity structure of the VF for the critical correlation functions, we must be careful 
about the relation between the two limits, g ^ and R — > 00, where i? is a typical IR cutoff. However, 
because of our proof for the finiteness of the VF for OPE coefficients, we can take the g ^ limit of the 
correlation functions before evaluating the integrals; the nested structure of the IR limits eliminates 
the kind of divergences which ordinarily appear as singular expressions. In section 5, we will verify this 
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point in the Ising model. 



4 UV Structure of VF for the Ising Model 



In the rest of the paper we concentrate on using VF for the calculation of various off-critical 
quantities for the two dimensional Ising model (IM). Since there are no dimensionless parameters in 
IM {i.e., the theory is super-renormalizable) the structure of UV divergences is relatively simpl^^. On 
the other hand, for spin fields the critical IM is not trivial and this makes the calculations of the finite 
parts and the structure of IR divergences of physical quantities very interesting and non-trivial. In this 
respect IM shares some of the complications of higher order calculations in other off-critical theories 
in higher dimensions and especially those in two dimensions near unitary minimal models. 

Singular OPE Subtractions 

For determining the necessary subtractions in the VF (with respect to m), we need the following 
singular OPE coefficients 

Omir)M0) = C^l\r;m)c^k{0) + O{\). (4.1) 

r 

From the general solution for OPE coefficients (l2.16| ), the singular part of Cm is given by 

C7(f)(r;m) = -G(r;m)H^\mr)G-\r;m), (4.2) 
r 

where the analyticity of OPE coefficients implies that Hrn {m) is a power series in m. This together 
with ( |2.35| ) and the fact that the singular part of the OPE coefficients must be at least as singular as 
^ lead to Hlr}{m)^ finite polynomial in m: 

H'nia\m) = "^"2 ' ^ dmHl^J'{()). (4.3) 

fc=0 

On the other hand, from the consistency relation ( |3.6D we have 

27rH^\m) = dm^{m) + Om{m) + mdmCm{m) + [cm{m),r{m)], (4.4) 

where we have used Om,g^(m) = 0. Inserting ( [4.4[ ) in ( [4. 2D , we can write the singular OPE coefficients 
in terms of the finite counterterms and anomalous dimensions. To calculate the subtractions we need 
to calculate the volume integral of the singular OPE coefficients. Note that the matrix G{r; m) satisfies 

^^For a non-trivial example, as far as UV divergences are concerned, consult ref. |^] in which the VF was used to 
study the perturbatively renormalizable ((/>^)4 theory. 
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r—G(r;m) = —G(r;m)T{mr). 
or 



(4.5) 



We can therefore write the singular OPE coefficients as total derivatives 



1 d 

2iTr dr 



G{r; m) [dm + rCmiiTir)^ G ^(r; m) + 7r/3imr^(lnr — -) 



(4.6) 



This means that we can calculate the integral in (^.41), and therefore the subtractions can be given in 
the following closed form 



Amie;m) = G{e;m) {dm + e Cmime)) G ^(e;m). 



(4.7) 



Hence, we have an explicit expression for the subtractions in terms of the matirx G{r; m) and the 
finite counterterms Cm{fn) only. Recall that the anomalous dimension matrix T{m) ( 2.35D fixes G{r;m) 
completely. Therefore, for the case of Ising model, we have completed the first step in the calculational 
scheme which was outlined at the end of section 3. 



Calculation of the Anomalous Dimension of the Conjugate Operator 



Let us consider the VF for the composite operator Om- From ( p. 21) we have 



dm (C^n 



lim 

e->0 



/ d\' {Om{r')Omm'in + ^mm'(e; m) 
Jr'>e 



(4.8) 



where from ( |4.7| ) and ( p. 35 ) we have the following expression for the subtraction 



Amm^ie; m) = -f3i ln(e) + Cmm^{rn). 



(4.9) 



Because of the relationship between the scaling limit of the IM and the theory of a free massive 
Majorana fermion, we have the advantage of having a nonperturbative answer for correlation functions 
involving only the energy density operator. The off-critical energy density operator is the most general 
scalar operator of scale dimension one, and in the fermionic representation it is a linear combination 
of i-ipijj and ml 



Cm = 7r~(^V'V') + ami. 
zvr 



(4.10) 



This fixes Om to be £" = ■^{iip'il>) at the fixed point. Now from the theory of a free massive Majorana 
fermion in two dimensions we have 1 23] 



{Om{r)Ommm = (ii^i^ir) ■ ii^mrm = (^^(mr) - K^{mr)) , 



(4.1i: 



22 



independent of a. The integral in the VF (4.11) can be done analytically and by comparing the result 
with the derivative of ( 2.27 ) with respect to m, we find the anomalous dimension of the conjugate 
operator to b^ 



(4.12) 



and the integration constant appearing in ( ^.27 ), ( 2.34| ), and the final expression for the free energy 
density ( p.28 ) to be related to the finite counterterm as follows 



1 1 

C = - exp[7 - 27rCmm J, 



(4.13) 



where 7 is the Euler constant. Therefore, upon fixing the finite counterterm Cmm^-, the integration 
constant is also fixed. 



Choice of the Parameters 

The precise definition of the conjugate operator Om depends on the choice of the parameters g\ 
and m which are by no means unique. We can reparametrize the space of renormalized trajectories by 
some general mixing of the old parameters up to two general restrictions |15]: The first is that the scale 
dimensions of the new parameters should be the same as the old ones, and be additively conserved 
under the RG flow. The second is that, in keeping with the locality of the theory, the reparametrization 
has to be analytic {i.e., it is given as a power series). So in our case we have only one mixing as follows 



~gi=gi + f (4.14) 
rh = m, (4.15) 

where we have kept the normalizations fixed. Under this reparametrization, the conjugate operator 
and the only relevant finite counterterm change simply as 



= Cm - ami, (4.16) 

Cmm — Cmm ~l~ Oi. (4.17) 

Hence, the mixing of Om with 1, or equivalently the change of parameter can be controlled by 
fixing the choice of Cmm^- Note that, as already discussed in section 2, the correlation functions do not 
depend on g-^. This and the fact that any change in Cffim can be absorbed in a change of g\, amount 
to the Cmm"*" independence of the correlation functions. 

If we consider off-critical correlation functions which involve only operators corresponding to spin 



and energy density operators, then for the possible subtractions needed in the VF, (f)a in (4.1) should 



^'^To find /3i we do not need to calculate the integral in VF. It is enough to compare (4.!:) with the corresponding 
subtraction coming from the following singular OPE coefficients |26[ 

^(s) 1 
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be only the spin operator a or the conjugate operator Om (which respectively correspond to a and £ 
at the fixed point). In this case, because of the dimensional restriction, the only relevant subtraction 
is the following which we have already seen in the VF for Om- 

Amm^{e;m) = ^ln{e) + Cmm^. (4.18) 
zvr 

This is the main reason that the necessary UV subtractions in VFs for IM are relatively simple, 
especially in higher orders of perturbation. Generally, this is a feature of super-renormalizable theories. 
Consider the case of the n— order VF. This means that we have to consider the insertion of volume 
integral of O^s, and the UV divergent subtractions needed when the inserted operators approach one 
another or other operators in the correlation function. But, the UV subtractions corresponding to 
two connected energy density operators are zero. This can be seen in the VF because of the obvious 
fact that: 1"^ = 0. In short; most, but not all, of the UV divergences in higher order perturbative 
calculations of IM are taken care of by the connectedness of the insertions. The rest of the subtractions 
require the calculation of the other elements of the anomalous dimension matrix T(g). /3i is the only 
element needed for the calculation of the short-distance expansion of the spin-spin correlation function 
to third order |2^. 



5 Perturbation Theory for the Ising Model 

In this section we implement the calculational scheme, outlined in section 3, for perturbative analysis 
of the spin-spin correlation function in the IM in the vicinity of the critical point. After studying the 
structure of the n— VF, we carry out the first two steps of our scheme for the n— order calculation 
of corrections to scaling. First, we study the structure of the necessary critical correlations. Then, we 
present a systematic method to calculate the necessary integrals over the critical correlation functions. 
The final step, which is the fixing of the off-critical operators and solving the resultant algebraic 
equations, is done up to third order. This section is a systematic extension of the explicit third order 
calculation done in reference p^], the details of which are presented in appendix B. We leave a systematic 
study of the final step for general n— order corrections to scaling, for a future work[p^]. 

5.1 n— Order Correction to Scaling 

Consider the short-distance expansion of the spin-spin correlation function. Due to the OPE algebra 
at the critical point, using the notation of ref.[^ we have 

(a(r)a(0))^ = a,^'W(r;m)(lW)™ + a/'W(r;m)(fW)™, (5.1) 

where l'^''^' and £^^^ are the off-critical deformations of the scalar operators in the families [1] and [£] 
(see appendix A). In the above summations, it is enough to consider only those operators which are 
not total derivatives, since the total derivatives have zero expectation value. 
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The OPE algebra of the IM at the fixed point has a Z2 symmetry, known as Kramers- Wannier 
duahty, which implies that the correlation functions involving operators in [1] and an odd number 
of operators in [£] are zero. As we have already argued at the end of section 3, this eliminates 
the contribution of the corresponding derivatives of the one point functions in the VF for the OPE 
coefficients. Hence, as a natural extension of the results of section 3, using 

a;;,(</)fc) =0 for n<Xk, (5.2) 

and keeping the subtlety in cancellation of IR divergences in mind, we obtain 



k even Xfj<k k odd x^<.k 

Note that there are a finite number of operators with Xf, < n, and therefore there are a finite number 
of IR subtractions. 

Using the short cut developed in section 3 we can write following finite and truncated n— order VF 



d:i,C^a\r)= lim ••• hm { {-irr"'-^Mn( 



r r 



-EE ^^'-^'-^'''^ 9r'=)a.i^W(r) • • • ^ 

k—l b perm {pi, - ■■ ,Pj^ } 

fc even x^<k from {l, --,n} 

+ E r^'-'^'^''^ d^r^-'^C^.'-^^'Hr) (5.4) 

fc — 1 b perm {pir'-jPfc} 

k odd a;j<fc from {1, ■■■,«} 

where the main term is 

Mn{Ui---Un)= f ••• / (funr''-^U£irUn)---£{rui)air)a{0)y, (5.5) 
and the IR subtraction integrals are 



'52fc,l;{6}(^l •••^2fc) 

lim • • • lim 

6^0 «fc^o 



d^ui ■■■ (fu2k r("M!.}+2fc)(^c(^n2fc) • • • ^^(rni)lW(O)) 

-'S2k^i-{b}{(.ir--(.2k,Ui,--- ,U2k) , (5.6) 

and similar expressions for the IR subtraction integrals ^^.^^i ^.^j^y The UV subtractions ^^^.^f^i 
and ^2^1 £.^i,y similar to equation ( 3.29| ), can be found in terms of Am{e), its derivatives, and fewer 
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volume integrals over correlation functions. The explicit forms of the VFs, up to third order, are given 
in appendix B. 

The VF ( |5.4| ) shows the nature of the communication between UV and IR divergences explicitly. 
The precise statement is that the perturbation for spin-spin correlation function does not have any 
explicit UV divergences, it involves IR divergences in which show up in the end as non-analytic 
dependence of the correlation function on the mass. To reveal this structure we need to handle the IR 
divergences by using the IR subtraction terms, involving '52^^.|^j and S2^^i ^.^f^y, which are naturally 

given by OPE. But S^^s have both IR divergences (which cancel the IR divergences of the main term) 
and UV divergences (cancelled by the necessary UV subtractions, <S^^s). 



5.2 Calculational Machinery 

Critical Correlation Functions 

The free massive Majorana fermion representation for the critical Ising model||^] makes the calcula- 
tion of the energy-density correlation functions trivial. Using the relation ( [4.10| ) at m = 0, calculation of 



correlation functions involving only the energy-density operators reduces to a product of Pfaffians[23|. 
Then correlation functions of operators in families [1] and [£] could be calculated by using the Vira- 
soro algebra. This takes care of all the correlation functions we need for IR subtractions. Note, as 
an important feature of these correlation functions which makes it possible to calculate their volume 
integrals, that they do not involve square roots. 

The expression of a{r) in terms of the fermion field, on the other hand, is non-local. This makes the 
calculation of the correlation functions involving the spin operator non-triviaf^. In two dimensions, 
due to conformal invariance, all correlation functions are calculable in principle Q. Moreover, in the 
case of the Ising model there are additional structures which help in principle to write systematically 
explicit recursive answers for the correlation functions of the type needed in the main term. Duplicating 
the IM gives a c = 1 theory which is amenable to bosonization, and hence the squares of the correlation 
functions could be calculated. This relationship was recognized much earlier as the equivalence of the 
Ashkin- Teller model at the end of the self dual line to two non-interacting Ising models] 24 1. Using 



the above relationships Di Francesco, Saleur, and Zuber found recursive relations for the correlation 
functions involving energy and spin operators [p^. For correlation functions appearing in the main term 



(5.5), in terms of a free bosonic field (p we have 



(r„) • • • £{n)a{r)a{0)) {a{r)a{0)) = -^icos 2(/)(r„) • • • cos 2</.(ri) cos ^{r) cos 0(0)) 

1 n— 1 

"2^ E {£irpJ---£{rp,^Mrym{£irp^)---£irp>{r)aiO)). (5.7) 

}c = 1 permutation of 
{Plr--,Pn} 
from {1, - ■ - ,n} 

See appendix A for the details and the calculation of correlation functions needed up to third order. 



From the above recursive relation and the neutrality condition ( A. 131 ) it can be shown inductively 



that the connected correlation functions have the following structure 

^*It is mainly due to this that even though the critical IM has a representation in terms of a free massless Majonara 
fermion, it is still a non-trivial fixed point. 
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{£{run)---£{rui)a{r)a{0)r = r^--^l[-—-. Hi ^ • ■ ■ , (5.8) 

^J-^ 27r\uk\\uk - 1| ^^^ti Fp - Uq\^ 



where F„(tti, • • • is a real analytic function without any square roots like \up — Uq\. Moreover, F„ 
has the following short-distance structure 

• • • 0(np - tig). (5.9) 

Hence, the connected correlation function does not have any unintegrable singularities. Therefore we 
have the following structure for Mn 



Ju-,>u^ 27r|Mi||Ui - 1| 7c„>„„ 2-K\Un\\Un - 1| ^\ti Pp " ""gl 

(5.10) 

The need to take the IR limits to infinity in ( |5.4| ), gives us a very important calculational advantage: 
usually it is much easier to calculate integrals for infinitely large cutoffs than for finite ones. This means 
we are looking for calculation of Mn and 5^^'s for 



K f/i < • • • < [/„. (5.11) 

Even after the above simplification, the calculation of A4„ involves elliptic functions, and beyond first 
order is extremely hard to do directly. In the following, we develop a systematic way to handle these 
integrations for any order. 

Symmetrization Transformation 

For each volume integration in 7W„ one copy of the integrand of 

M,{Ur) = f ^'^^ (5.12) 

Jui>ui 47r|ui||ui - 1| 

appears. According to the above analysis this is the only source of square roots. Therefore, we can 
pinpoint the root of the calculational difficulty in the very first order integral A4i. The key point 
in overcoming this difficulty is to use the symmetry properties of the integral. This can be seen by 
comparing the integral in (5.5) with the ones in the IR subtraction (^.6[); in the IR subtraction terms 



we have both IR and UV divergences. More precisely, the structure of the IR subtractions is invariant 
under the simultaneous inversion of all variables of the integrals. In the case of the main term only the 
A^i-like factors do not respect this symmetry. The absence of UV divergences in the main term ( ^.5| ) 
is directly related to separation of the singular points in the Ali-like part of integrand, one at zero and 
the other at 1. Therefore it is enough to look for a transformation which symmetrizes factors of this 
type. 



27 



This desired transformation changes the complex variable r j with some linear function of Xj + l/xj , 
where Xi is the new complex variable. Consider the conformal transformation from Xi to 



(5.13) 



r- plane 




/ 



lr'-r"l 
4R 



X— plane 




1 4R 

lr'-r"l 



Figure 1: Symmetrization Transformation 



which is called the symmetrization transformation! 26]. It merges the two integrable singularities in r[ 
and r'l to one unintegrable singularity at zero and enforces the symmetry between UV and IR. We 
only need to calculate the integrals for infinitely large cutoffs as it is given in ( |5.11 ). This is crucial in 
having a simple transformation for the IR cutoffs under the symmetrization (See figure 1). Moreover, 
using the inversion symmetry which has been produced by symmetrization, we can restrict the domain 
of integration to the region with radius larger than 1. So we have the following transformation for a 
typical integral under symmetrization 



(fxj 



(5.14) 



By applying the symmetrization on eg. ( p. 10 ) we get the following expression for Ain 



Mn{Uu---,Un) 



cfxi 



4,ui>xi 27r|a;i| 

n 

n 



d Xyi 



Xn 



XpXq 



12 H 



(xi + 1)^ 
4x1 



{Xn + 1)^ 
4Xr,, 



(5.15) 



where we ignore the irrelevant higher order corrections because of the condition ( |5.11| ). The sym- 
metrized integrals are much simpler than the integrals we started with. Their integrands are rational 
functions, and they can be calculated by standard methods. For explicit calculation to the third order 
see appendix B. 
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6 Spin-Spin-Energy Density Correlation Function OfF-criticality 



In this section we demonstrate the apphcabiUty of our method for the calculation of more general 
correlation functions. As a typical example we calculate the short-distance expansion of the spin-spin- 



energy density correlation function to first order. Our starting point is equation (2.11) for n = 3 



{£{r')a{r)a{{)))m = C^^Jir' , r; m)(</>b)^, (6.1) 

where in the sum over h it is sufficient to include scalar composite operators, which are not total 
derivatives. Because of the operator algebra of the critical IM, G {[1], [<S]}. 

Zeroth Order: The Scaling Behavior 



At m = 0, due to equations ( pOj) and ( |Al^ ), the OPE (|3) b ecomes 



3 

a.n^'(r,rO= . , „ , . (6.2) 
47]- r' — r r' 



First Order Corrections to Scaling 



The first order variational formula for general OPE coefficients (3.14) can also be truncated. The 



resulting variational formula can be brought into a more convenient form in which the main term is fully 
connected. This is motivated by the results of section 4; the necessary UV subtractions are separated 
out of the main term, and therefore our integration methods are applicable. We can make the main 
term fully connected by adding and subtracting an appropriate term. Therefore, using the short-hand 
notation for truncated VFs, we have 

dmC.„ra^(ry) = 9^ (r')a(r)a(0))^ + ((a(r)(7(0)) - C^am'^{r,r'))dm{£), (6.3) 
which according to the calculational scheme of section 3 is equivalent to 



a™a<x™^r,r')= lim (r-^ M (-;—)- (^.^'"(r, r') - r"!) (5i,£( — ) + c^^i) 1, (6.4) 



where iSi^^^ is defined in (|B.5|) and 



Mi{u'-U")= / (Pu" r-i {£{ru")£{ru')(j{r)cj{{))) 

Ju">u" 

(fiu" 



1 



47r u' ki' — 1 Jiu">u" 47r m" u — 1 



u' + u" - 2u'u"P - A\u'\\u' - l\\u"\\u" - 11 



\u' — n"P 



.(6.5) 



After the symmetrization transformation for infinitely large U" we obtain 
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4x 



41x1' 



7r|x2 - Ip Jiu">y>i 47r|yp 



Re 



{x - y){l - xy) 
{x - y){l - xy) 



(6.6) 



This integral can be calculated with standard techniques. With our previous choice of finite counterterm 
Cmm^, the final results are 



27rr 4 



|9».' ».|2 

'^aam V i ' / — 1 i 

4|?^'|[t'' — 7*17^4 

(1 — r P \ 
,,, , ' , + 1 ) Inr - ln(16|r'||r' - r\) + 2- 
4 r \r — r ^ 



x' - 11 



In Ixl 



(6.7) 
(6.8) 



where x = j:{2r' — r ± 2\Jr'(r' — r)). Now inserting back into ( |6.lD we get 



3 

7" 4 



47r r' — r r' 



1 + J ^ ln( mr e^/8 

2r 



2|r' — r| |r'|m J /16|r' — r||r'| \ mr / 
r V r / 4 V 



In l^l + 0(mV2) .(6.9) 



The two different short-distance limits of this result can be easily checked to be consistent with those 
obtained from two-point functions already calculated. This is a typical example of the applicability of 
our calculational scheme for non-trivial higher point functions. 



7 Conclusion and Future Directions 



Three things were done in this paper: First we reviewed the new formalism for the study of 
the short-distance structure of renormalizable field theories. This formalism provided us with the 
first systematic scheme for perturbative calculations near a non-trivial fixed point. This review was 
done both in general and for the special case of the two dimensional massive Ising model at zero 
external magnetic field. In particular, we presented a complete renormalization group analysis of all 
the one-point functions of the theory which contain all non-analytic dependence on the mass of any 
correlation function in the theory. Second, we provided the details of the proof of finiteness of the 
variational formula for operator product coefficients to all orders. This amounted to the consistency of 
the variational formula with the analyticity of the OPE coefficients with respect to parameters of the 
theory. As a by-product, we developed a systematic scheme for perturbative calculations of the OPE 
coefficients in the vicinity of any fixed point. Our first main result was the truncated and manifestly IR 
finite VF (3.24), in which there were only a finite number of terms. This VF can be extended to higher 
orders readily. Third, we implemented this scheme in the analysis of the short-distance expansion of 
the Ising model correlation functions which only involve two spin operators. The truncated VF for 
the case of spin-spin correlation function to n— order was given by equation (0). In this case we 
developed a systematic method for calculating the necessary integrals. Our second main result was 
the symmetrization transformation (5.13) which is at the heart of our method and transforms these 
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extremely hard integrals into a managable form. In particular we calculated the spin-spin and spin-spin- 



energy density correlation functions to third and first orders, equations ( B.16 ) and ( |6.9| ) respectively. 
Our result for the spin-spin correlation function is in complete agreement with the scaling limit of the 
well-known exact calculation on the lattice. 

This calculational scheme is formulated in coordinate space and its main advantage over con- 
ventional schemes is its applicability to perturbative analysis of renormalizable field theories near a 
non-trivial fixed point. The fundamental strength of this formalism relies on a natural separation of 
infrared and ultraviolet divergences by incorporating OPE from the outset. An important subtlety 
to bear in mind is that despite the non-perturbative validity of this formalism, its applicability to 
perturbative calculations is much more intricate for theories with some ill-behaved conjugate operators 
at the fixed point However, our scheme is well-suited for a large class of two dimensional field 
theories, defined as some perturbation of unitary minimal models ^, One of the difficult steps in 
the analysis of these models is to calculate non-trivial finite integrals of the critical correlation functions 
in each theory. We hope that the possibility of calculating similar integrals in the Ising model sheds 
some light on the structure of integrals in the more complicated minimal models. Clearly, there is 
much interesting work to be done. 

Among the possible future directions for extending this formalism, three stand out. The first is the 
renormalization group analysis and perturbative calculation of the correlation functions in the Ising 
model with non-zero external magnetic field. The second is to understand the structure of the flow 
between tri-critical and critical Ising models. The third is the completion of the final step of our 
program for the Ising model, as outlined in section 3. It would be extremely instructive to explicitly 
see that the validity of our scheme to any order fixes all the anomalous dimensions and the finite 
counterterms. It should also clarify the role of integr ability of the massive Ising model in organizing 
our field-theoretic calculations effectively [p^j. 
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A The Critical Ising Model 



Operator Content 

At the critical point the Ising model acquires conformal symmetry, and we therefore analyze the 
operator content of the model using the techniques of conformal field theory (Consult [0] for more 
details.). The local operators of the critical Ising model fall into three classes: class [1], [£], and [a]. 
Each class is denoted by their respective primary operators. We call an operator primary if 



Ln(p = Ln(p = for n > 0, (A.l) 

where Ln, Ln{n G Z) are left- and right-moving Virasoro operators. In the class of the primary operator 
<j) a general operator is given in the following form 
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^W{n} ^ . . . L_^^^l_^_ . . . L_^_^. (A.2) 
To find the linear relation among various operators we need the Virasoro algebra: 

[Lm, Ln] = (m - n)Ljn+n + -^^{im? - l)5m+nfl, 

[Lm, Ln] = (m - n)Lm+n + ■^m{im? - l)5m+nfl- (A.3) 

The scale dimension x and the spin h of the local operator 0{"j"t"}j are given by: 
X = (A + niH +nAr) + (A + njH I-^tv)' 

h = (A + ni + --- + n7v)-(A + ni + --- + n^), (A.4) 

where (A, A) are the eigenvalues of {Lq,Lq) for 0; for 1, and a they are given by (0,0), (-, -), and 

(— , — ) respectively. Therefore, for scalar operators a lighter notation than the one in ( |A.2D can be 
adopted; = </.^"H«}. 

The Ward identity associated with SL{2, C) invariance of conformal field theories implies 

L-i =d, = {d, - idy)/2. (A.5) 



Therefore, using (A.2) and (A.3) we find that, up to scale dimension 7, we have the following scalar 
operators 

Oo = l, 0„, = £, T = 1^2} = L_2-^-2l, (A.6) 

in the families [1] and [£]. 

The non-trivial operator product expansions between members of these conformal classes are given 

by: 

[a].[a] = [l] + [£], [a]-[£] = [a], [^] • [^] = [1]. (A.7) 

The off-critical Ising model is obtained by adding local scalar operators to the critical hamiltonian 
density. Therefore the most general hamiltonian for the off-critical Ising model is 

^J■[£], fJ-ia]} = J J J (A.8) 

where the term /xjjjOjj] denotes a linear combination of all scalar operators in class [i]. By virtue of 
(A.7) and the fact that RG transformations only depend on local physics we can set = 0. Note 



that, since the OPE of two operators in class [a] contains no operators in that class, it make sense to 
speak of off-critical correlation functions with an even number of spin operators. 
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Correlation Functions 

We start by fixing the normalization of our operators at the fixed point 



i^i^nO)) = (A.9) 
(^(0-(0)) = (A.10) 

The convention that correlation functions have no subscript means that we are calculating fixed point 
quantities. 

We now use the results of to calculate all the necessary correlation functions of the critical Ising 
model. All the correlation functions are expressed in terms of the correlators of a free bosonic field (j) 
whose two point function is 

{(j){r)(j){r')) = -]^hi\r -r'l (A.ll) 

which implies 

mn 

Jm4>(r) in<t>{r') 



and the neutrality condition 



^■Q (.imiHn)^ ^0 iff ^ mi = 0. (A.13) 

i=l 1=1 

The non-trivial three point function is given by (See equation (2.27) of ref. 

3 

{£{r')a{r)a{0))o = -(a(r)a{0)) ■ (cos </>(/) cos (/)(r) cos (?!)(0)) = — , ,, ,, . (A.14) 
vr 47r r' — r r' 



By letting r — > and using ( [A.9| ) and ( [A. 10 ) we get: 



3 



a<x'"(r;0) =7rr4 (A.15) 



The four point functions of interest are {££aa) and {££££). Similar to ( A.14 ) we start with (See 
equation (2.28) of ref. 



{£{r")£ir')air)aiO)) ■ (a(r)a(O)) + (f (r")a(r)cj(0)) • (^(r')cj(r)a(0)) 

= (cos 20(r") cos 20(r') cos (/)(r) cos (/>(0)) (A.16) 



Note that our particular choice of normalization in (|A.9|) modifies the equations from ref. l25| by factors of 2tv. 
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and using (A.14|), (A. 10), and (A. 13) we have: 



{£ir")£{r')a{r)aiO)) 



I'P'p^^ — |— ipip^ 2 7^ " T" ' I ^ 



(A.17) 



from ref. [25| we also have: 



E 



cyclic pcrmutatio: 



|ri|^|r2 



kiPk2 



(A.18) 



The last correlation function we need is {£££aa). As pointed out in |25|, a slight extension of 
( A.14| ) and ( [A.16| ) allows for calculation of correlation function of any even number of spin operators 
and an arbitrary number of energy density operators, we therefore have 



{£{r3)£{r2)£in)air)am ■ (a(r)a(O)) + (£:(r3)^(r2)a(r)a(0)) • (£:(ri)a(r)a(0)) 
+ {£{rs)£{n)air)a{0)) • {£{r2)air)a{0)) + {£ir2)£{n)air)a{0)) ■ {£{r3)a{r)a{0)) 

= 8(cos 2(/)(r3) cos 2(j){r2) cos 2(f){ri) cos (/)(r) cos 0(0)) 



(A.19) 



from which, using ( |A.17D , (|A.14|) , (|A.iq) , and ( |A.13D , we get 



(f(r3)f(r2)£:(ri)a(r)a(0)) 
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r 4 



167r'^|ri||ri 

X E 



cyclic permutation 
°t {''3''"2'''l) 



rp|r3p 



|r3 — riP|r3 - r2p 



|rri + rr2 — 2rir2| 
4|ri — 



.(A.20) 



More general correlation functions can be calculated similarly. 



B Calculation of Spin-Spin Correlation Function to Third Order 



In this appendix we present the details of the calculations of ref. [26|. This is an application of the 
analysis of the n— corrections to scaling for the spin-spin correlation function (section 5). 



Zero-th Order: The Scaling Behavior 



As a consequence of ( 3.20| ) and ( [A.IOP we have 



a.i(r) 



,1/4 ■ 



(B.l) 
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Corrections to Scaling 

For n = 1, 2, 3, the equation ( ^.4D reduces to 




lim lim lim {r ^ -M-zi — ~) — ~i — ~1 

\ f J' f J 



_Ri— >oo i?,2— »oo Rs-^oo 



' ^ ^ ^ J 

cyclic perm. 



cyclic perm 



where the expression for Aln, (^^ = 1)2,3) is given by (^.51) and the subtractions are as follows 



= lim f / J^p—^—\^[^ 



(B.5) 
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The explicit expressions for necessary correlation functions are given in appendix A. For the calculations 



of the main terms we use the symmetrization transformation ( 5.13 ). The resulting integrals are 



M2{Ui,U2) 



4Ui >ii >1 



MiiUi) 
Att\xi |2 
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47r|x2|' 



Re 
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(B.7) 
(B.8) 



>t3(C/i,f/2,C/3) 
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4Ui>xi 27r|xi|2 J UJ2>X2 2vr|x2|2 J4u^>x^ 27r|x3|2 |x„ — x„|'^|l — x„x 



P<1 



X] jl^l ~ ^2^1 - 3:iX2|'^|x3 - l]'' - |xi - X3p|x2 - a;3p|l - X1X3PII - X2X3P 

Ip + 2i?e((xi - X2f{l - xiX2)^)) |. (B.9) 



cyclic 
perm 



These integrals do not have any square roots and they are easily calculable. These allow for the 
calculation of corrections to the scaling behavior of the OPE coefficients. The results are 
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11 
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■ In r. 



(B.IO) 
(B.ll) 
(B.12) 

(B.13) 
(B.14) 
(B.15) 



The results in equations (B.ll) and ( B.15| ) are due to the particular choice of the finite counterterm 
Cmm"*"- Here Cmm^ = ^ ln2. Consequently from ( 4.13| ) we fix the integration constant c = where 7 is 
the Euler's constant. Recall that the answer for the correlation functions is independent of the choice 
of the finite counterterms (See section 4). With this information at hand, we arrive at the following 
expression for the short-distance expansion of the spin-spin correlation function 



1 



T 4 



1 



a{r)a{0))m = ^{l + ^t H\t\ey8) + + ln(|t|eV8) + 0{t^ In' \t\)), 



16 32 



(B.16) 



where t = mr. This is in complete agreement with the results of the scaling limit of the lattice 
calculation P]. 
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